Preliminary results on the homogenization of thin 
piezoelectric perforated shells* 



Houari Mechkour^^ 



Abstract 

We consider a composite piezoelectric material whose reference con- 
figuration is a thin shell with fixed thickness. In this work, we give a new 
approach based on the periodic unfolding method to justify the modelling 
of a thin piezoelectric perforated shells and we establish the limit consti- 
tutive law by letting the size of holes is supposed to go to zero. This allows 
to use the homogenization technique to derive the limitting equations and 
the homogenizaed coefficients are explicity described. 
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1 Introduction 

The shells is a three dimensional continous medium, where its thickness is small 
compared to other dimensions. Its geometry is characterized by two small pa- 
rameters : the tickness of the shells and the size of perforations. The behavior 
of piezoelectric shells when the tickness goes to zero has been studied by Haenel 

In this work, we consider a periodically perforated piezoelectric shells. For 
two dimensional limit equations obtained by Haenel [S], we study the behavior 
for the elastic displacement and electric potential as the size of perforations 
becomes smaller and smaller. We are concerned with two independent problems 
: the membrane and bending problems and we give the convergence results 
based on new periodic unfolding method, recently introduced by Cioranescu, 
Damlamian and Griso 

We find explicitly the overall homogenized tensors and we study their prop- 
erties in order to give a corrector results associated for membrane and bending 
problems. 



*This preprint is translation of Chapter 3 of PhD thesis [9] (Full text: 
[iittp: /7tel.archives-ouvertes.fr/tcl-00008496/fr/ ) 

^Centre de Mathematiques Appliquees (UMR 7641) Ecole Polytechnique, 91128 Palaiseau, 
France, (mechkour@cmap.polytechnique.fr). 



1 



2 The periodic unfolding method 



The Periodic Unfolding Method is a novel approach to periodic homogenization 
problems that applies as well to problems with holes and truss-like structures or 
in linearized elasticity. The periodic unfoding method is equivalent to two-scale 
convergence method but is both simpler and more efficient. 

First we briefly introduce this method which was developped in [5] . 

Let ri C be a domain with a smooth boundary 9il and let Y = [0, 1]'^ be 
the reference cell. We consider that S is an open smooth boundary subset of Y 
such that S <ZY and set Y* —Y \ S. S plays of the reference hole while Y* is 
the part of Y occupied by the material. We also set S'^ = e{S -I- fc) n f2, A: e Z'^. 

The perforated domain Vl^ is defined as the set fi^ — il\S'^. We assume that 
17^ is connected and that the holes do not intersect the boundary dfl. 

For z g IR'^, we denote by [z]y' the unique integer combination such that 
z - [z]y' € Y* and set {z}y- = z - [z]y. e Y*. 
For any x G M"^ and e > we have 



X ~ e 



+ {-} 



Define : L^{n^) L^{n'' x Y*) with 

'X 



T^{w){x,y) =w{e 



for all (x, y)en' X Y* 



Obviously, for any e _L^(r2^) we have 

T'iv + w) =r''{v)+T'{w) 



(2.1) 
(2.2) 

For our purpose, all functions defined in L^{Q,'^) are extended by zero outside 

Proposition 1. (Properties of T^) (see [2]) 
(a) For all w G L^(ri'^) we have 

1 



wdx 



Y* 



{w)dxdy; 



(2.3) 



(b) For any w G L'^{J¥) we have 

T^iw) w strongly in L^{Q xY*); 

(c) Ifiw^) C L^{n^), then 

^ w weakly in L'^{ri) T^{w^) w weakly in L'^{fl x Y*); 

1 



(2.4) 



T^w^) w weakly m L\VLxY*) 



w 



Y* 



wdy weakly in L (Q). 
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Proposition 2. (see Let (w^) C L^(ri^) be a bounded sequence. Th 



T^{w^) w weakly inL^{fl x Y*) 



w"^ two-scales converges to w. 



Theorem 2.1. (see [2]) Let (w^) C H^in") he a bounded sequence that weakly 
converges to w in H^{n). Then there exists £ {ft; H^^^lY)) such that 

T^w^) w weakly m L^{n x Y*), (2.5) 

T^V^w^) V^w + Vyw^ weakly in L^{VL x Y*). (2.6) 

Theorem 2.2. Let {u'^)e be a bounded sequence in H^{fl^). Then there exists 
u e H^{n) and € L'^iVl; H^^^{Y)/R) such that 



T^{u'^) u weakly in 

T^(S7xU^) ^ Va;it weakly in 

T^Vlu"") Vlu + Vlu^ weakly in 



L^{n X r*), 

L^{n X Y*), 
L^n X Y*). 



(2.7) 



Proof. Since (u^) is bounded in H^{fl'^), it follows that it weakly converges to 
some u e 

ij2(17). According to TheoremO there exists G L^(rJ; H^^^{Y*)) 

such that 

T^iu") u weakly in L^{n x Y*), (2.8) 

and 

T''{V^u^) ^V^u + VyU^ weakly in 2.2(17 X y*). (2.9) 

Moreover, by the boundedness of (u^) in H^{il^), it follows that T^CV^u"^) is 
bounded in L^{^ x Y*). Hence, there exists g € L'^{fl x Y*) such that 



weakly in L^(nxY*). 



Let V G C^c,{Y*)). Then we have 



dxjdxi 



da; 



9-0 1 
dxi £ dyi 



Using the unfolding operator and 



in the above relation we have 



V dxidx-i 



il){x,y) dxdy 



dxi 



dip ^ dip 
dxi dyi 



Passing to the limit in (|2.12p and using 



(Pl^ we get 



(2.10) 
(2.11) 

(z, y) dxdy. 
(2.12) 



= 



This yields that — h — — does not depend on y. Since is y*— periodic in 



dxj dyj 

the second variable we conclude that u^(x, y) — u^{x) and by 
r^(V^u^) ^ V^u weakly in L'^{Q. x Y*). 



we deduce 
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Let now ip £ C^^^.{Y*)) with Vyi){x, y) = 0. From (piT|l we obtain 
, / x\ , f du^ dtp 



dxidxj 



da; 



dxj dxj 



dx. 



Using again the unfolding operator we have 



dxidxj 



Passing to the hmit we get 

gij{x,y)ip{x,y)dxdy 



r dth 
^ix,y)dxdy = - J i^—j ^ix,y)dxdy. (2.13) 



OxY* 



The above relations lead us to 



OxY* 



dxidxj 



{x,y)ip{x,y)dxdy. 



Qt] {x, y) 



dxidxj 



ip(x, y)dxdy = 0, 



for aU V' G C^^^{Y*)) with Vy'4}{x, y) = 0. 

Then, there exists u £ [L'^{n;H^^^.{Y*)/R)]'^ such that 



(2.14) 



Taking into account the symmetry of the left side member in the above equality, 
we may conclude that there exists £ L^{^', Hp^^.{Y*) /M.) such that u — \7yU^. 
Now (I2.14P becomes 

The proof of Theorem 12.21 is now complete. □ 
We introduce now the averaging operator 

W : L^{n' X Y*) ^ L^{n'), 





x' 











+ ez, 
Y* ■> Y' 



ill 



dz. 



Proposition 3. (see [2]) 

(i) (j) strongly in L^in), for all (j) £ L'^i^); 

(ii) WiT'icf))) ^(P, for all(l)£ L^{n^); 

(iii) T^u^m^^y) = 



1 



Y* 



+ £z, |-| ^ dz, for all $ £ 
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(iv) [ W{^){x)dx = [ ^{x, y)dxdy, for all $ e L'^i^" x F*); 

(v) ^ *(^' 2/)^^' V a" $ e X y*). 

Theorem 2.3. (see [2]) Lei {(tf)e C T/ie following weak convergences 

are equivalent 

(i) T^(0^) ^ $ in X r*); 

(ii) (j)^ ^ (j) in L'^{n). 

A similar equivalence holds for strong convergences. 

3 The membrane problem 
3.1 The convergence results 

Latin indices take their values in the set {1,2,3} and Greek indices take their 
values in {1, 2}. The sumation convention is also used. Boldface letters represent 
vector-valued functions. 

Let cj C be an open bounded and connected set with a Lipschitz- 
continuous boundary. Let ip : uj ^ 'S? he a, C^(w) one-to-one function such 
that the vectors ai ~ dip/dxi and a2 = dip/dx2 are linearly independent at all 
points in uj. Denote 

^1 ^ ^2 , , , , s 

aa ~ r and a — detlao, • slr). 

|ai X aal 

We consider the membrane problem 

(3.15) 

for all (v, V') G Y{lu^) x W^(w^), where 

V(w^) = e H\lu') X H\uj^) X L^{uj^); t;„ = on 7^'^}, with meas > 0, 
W{uj') ^ {ijj e H^Lu"); V = muj'\uJo,TP^O on 7^}. 

c^(u, v) = c^'S^'^-^V^ 7a/3(u)7A^(v), 
e^(v, ^) = e^'^^'^V^ dx^Pj^piv), 

Sa/3(v) = i (daVp + dpVa) , 

=F\xi,X2)^ J J\xi,X2,t)dt, where f = (f ) e L2(c^ x [-1, 1]), 
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q=(gOeL2(r+ur_). 

Taking into account the ellipticity and the symmetry of the tensors {c"^^^) 
and {d"^) we can apply the Lax-Milgram Theorem in order to get the existence 
of a unique displacement u"^ and of a unique electric potential (p"^ that verifies 
the variational equation p.lSp . 

Theorem 3.1. The sequences (T'^(u^))^ and {T^{lp^))^ weakly converge to 
u e and ip G Wiuj) respectively, which are the unique solutions of the 

homogenized problem 

[c{u,w)+e{v,ip)]dxdy = \Y*\a [ F'v.dx + \Y*\a I q'v^dV, 

-e{\i,v) + d['ip,Lp)\dxdy =0, 



(3.16) 



for all (y,i]j) G V(a;) x W{lu), where \Y*\a — / \/ady and 



c(u,v) =c"'5^'^V^7a/3(u)7AM(v) 
e(v,^/>) =e^"'3v^9A^7a/3(v), 



pcrV // ' r ' ^ pcrV 
-a(3re ^ 

-cra/3 _ 



for all u,v e Hi,,(y*)/M, ^> e i/i,,(r* 



The local functions (w'^^,("^^) anrf {x"^ ,rf) verifies the local problems 
(c,(I]^''+w^«,v) + e,(v,r«))dy =0, 

(-6,(2^" + + ^)) = o, 

{cy{'L'',w)+ey{w,y„+-q''))dy =0, 

(-ey(z'', V) +?7'^,'0))rfy =0, 

for all {v,iP) e Hi„(r*)/R x Hl^^{Y*)/R, where S"'^ ^ y^e^ + y^Ca and 

c^(w,v) = c"'3-«VHs„^,^(w)sA^,j,(v), (3.20) 
e„(v,^) = e^'^'^V^SapMdxM, (3.21) 

d„(l/;,(/.) = d"^V^9a,yV^9A,y0. (3.22) 



(3.17) 



(3.18) 



(3.19) 
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Proof. To obtain an a priori estimate, we take v = u*^ and — ip'^ 'm p.l5p . 
By Korn and Poincare's inequalities for perforated domains, it foUows 

l|u''||v(t^-) + W\\w(uj-) < C, 

where C is a positive constant that depends only on uj (but not on e). 

By Theorem l2.11 there exist (u, ip) G V(a;) x W{uj) and two fields of corectors 
ui = (ul^ul) e L2(tj;i/i^^(r)/M), ^1 e L^{uj-HI^^{Y*)) such that, up to a 
sequence we have 

T'^(u^)^u weakly in L2(w X (-3 23) 

T^(V:r(u^) V^(u) + VyU^ weakly in L2(tj X y*;]R2)^ (3.24) 
T^V^{tp^) ~^V^^{(p) + Vytp^ weakly in L2(w X r*;M2). (3.25) 
The linerity of implies 

r"(7„;3(u")) ^7a/3,.(u) + Sa,3,y(ui) Weakly in (o; X ; M^) (3.26) 

and 

T^dap') da.xV + da.^yip^ weakly in L^{uj x Y*-M?). (3.27) 
We chose as a test function in (|3.15p 

v^(x) =vi(a:)+ev2 (x, ^) , vi e 2?(c^); V2 € 2?(c^; Cp°^,(r*)), 

V^^x) = V;i(x) + (x, ^) , ^1 e I?(a;); V2 e I?(^; Cp°S,(r*)). 
It follows that 

v''(j;) ^ vi(a;) strongly in L2(w), (3.28) 
S/x^^{x) VxVi{x) + V2 {x,y) weakly in l/{uj x Y*), 
T"{lp<y{^")) ^ IpM) + ep^.y(V') weakly in h'iuo x F*), (3.29) 
T%daip^) ^ da.^ipi + aa,;;^'^ weakly in L^{lj x F*). (3.30) 

Using the linearity of and the fact that [ v — r / '^^{v) for all 

Juj K *l JjuxY" 

V £ L'^{uj), we can pass to the limit in the variational form p.l5[) . We get 

/ [c°'^''^y/a {jxf,,x{u.) + SA^,y(u^) + e^^^y/a {dxif + dx.yP^)] (7q/3.x(vi) + Saf3.y{^2) dxdy 

^\Y*\a f F'vudx + \Y*\a f q'vudT+ f /i"''V^[(7a/3,.(vi) + s„^,y(v2)]da:dy, 
J Juj Jr_ur+ JujxY* 




I'^y/a {da.xi!! + da,y^2) dxdy 



(3.31) 



7 



Letting V2 = 0, '!/'2 = in the above equalities, by density we have 

- \Y*\a I F'v,dx + \Y*\a f q'v.dT + f /i"^^/^7,^,,(v)dxd2/ 
[ [-e"^'^V^ (7am,.(u) + SAM,y (u') + (9a <^ + dx,y<t>^)] d^^^ijjdxdy 

JuxY* 

1°' \fada^xi^dxdy , 



(3.32) 



for all (v, G V(a;) x W^(ijj). By linearity, we take 

M^{x, y) = 7re(u(a;, y))w^^(j/) + a„(^(a;)z''(t/) + y), (3.33) 

^i(a;, y) = 7re(u(a:, y))r'(y) + dMx)v''iy) + v). (3.34) 
where w"^^ , CJ^ , z'^ , 77'^ are Y* periodic in y and 

c^f'^^xMx, y)) + e'"''9A,,^(a;, y) = h^^ix, y), (3.35) 

~ e"^^sx^Jq{x, y)) + d'^'^dx^y^x, y) = y). (3.36) 

Replacing p.33p and p.34p in (|3.32p and taking into account the properties of 
q{x,y) and ^(x,y) we obtain 

/ [(f^^'jre{u)+r''^d„^] jc,p{v)dxdy = \Y*\a I F'v^dx + \Y*\a I q'vdT 

/ r —^OltO , , — Q(T - 

/ \-S 7re(u)+c( d\ip\da%jdxdy = 

Ju)xY' 

(3.37) 

for all (v, ?/;) G V(ti>) x W(uS), where 

c"^"' = ^ [c"''^'^A/^((5I^ + SA^,j;(w^«)) + e^"''y^aA,,r']dy,(3.38) 
e'^"'' - ^ {c^f'^^^rasy.^{7F)^e^^f'^a{b'i^dx,yTf)\dy, (3.39) 
7"^' = [e"^'^V^(<5I^ + SA^,,(w-«))-d"^y^9A,,r]dy, (3.40) 

[-e"^''V^SA^,j,(z'^) + rf"^V^((5J + dy,,yTf)\ dy. (3.41) 



We will prove that e'""'^ = f . In that sense we first determine the local 
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problems verified by the functions w""^, C^^, z"^, 77"^. From ([2351) -([Sill]) we have 
c"^"' = J [c"^'^^V^SAM,y(S"' + w"'') + e^"''V^aA,;;r']rfy, (3.42) 
r^'' = [c"^^^V^SA^(^") + e^"^V^9A.,(y.+77")]rfy, (3.43) 

T^' = ^Je"^^^^SAM.y((S^'+w^«))-d"^^/^aA^,rVy' (3-44) 
= [-e"^''V^SAM.y(^") + rf"^V^aA.j;(y.+77")]d2/, (3.45) 

We let vi = and ipi =Q in (|3.3ip . By density we deduce 

[c'^f^^^/^ {jx^A^) + SAM,y(u')) + e^"^V^ {dxv + dx,yip^)] sc.p.,y{^)dxdy 

xY' 

= / h^^y/a Sa[3,y{'y)dxdy, 
I [-e"^^V^ {iXt^A^) + SAM.y(ui)) + d"^^ {dxip + dx.yif')] do^^yHxdy 
I" y/ada,yipdxdy , 



xY' 



for all (v,V') e L2(c.;Hi„.(y*)/R) x L2(^. 
Using again p.33p and (|3.34p we deduce 

(c,(S-'' + w-^v)+e,(v,r^))dy = 0, 
(-e,(E^«+w-«,^)+d,(r«,V^))dy =0, 



(cy(z'',v) +ey(v,?/^ +77"^))^^ =0, 
(-ey(z",V^)+dj,(y, +77^,1^)) =0, 



(3.46) 



(3.47) 



(3.48) 



y. 



for aU (v, V) € Hi„.(r*)/M x iJi„.(y*)/R, where Cj,, Cj,, are given by ((3:201) - 
Next, we prove that e'^^" = 7'^^". By taking (v, t/j) = (z"^,??"^) in ([233), we 

get 



(c,(w-^z'^) + e,(z^r'))dy =-/ c-^^'^'^y/^sxA^n. 

Jy- 

{-eyiw-'>,r,n+dy{C',w))dy = / e^^" V^dx,yV^ dy , 



lY* JY* 

and so, by p.39p it follows 



(e'^-'' - c,(w-^ z-)) - ey{v,^\rf) - e,(z^ C^) + dy{^%C')) dy. 

(3.49) 



In the same manner, by taking (v, ?/;) = {z'^, r]") as a test function in (|3.48p . we 
get 

/ (c,(w^^z^) + e,(w^^77"))rf2/ =-/ e-^''v^SAM(w^''), 

V Jy JY' 

and by (|3.40p it follows 

(3.50) 

Now, by and we can conclude that e"^^^ = J'^^". 

To prove that the homogenized problem is well posed, it remains only to 
show the symmetry and coercivity of (c"^'^^) and (d ). 

(i) The ellipticity and symmetry of the tensor (c"''^^). 

Let us prove first the symmetry. From (|3.38[) we deduce 

It remains to show that c"^"^^ = c^^°'^ . From p.42p we have 

^P^' ^ J [c^,(S"«+w^^S"'')+eJ,(S"'5,C'*)] dy. (3.51) 

But 

(3.52) 

On the other hand, by taking (v, ^p) ~ {w"'^ , C,"^) in (|3.47[) . we get 

/ e,(w"^CVy c,(S-'' + w-^w"'3)dy, 

Jy- Jy' (3.53) 

lY* JY' 

Replacing now p.52[) - (|3.53p in (|3.51|) we deduce 

^a/3re^y Cj,(I]^^+w^^I]"''+w"'^)dy + y" dy{C'^,C^)dy. 
From the above relations we can easily conclude the symmetry of the tensor 

Let us prove the coercivity. Let (Xa/3) be a symmetric tensor (i.e. Xap = 
Xpa)- First we note that by p.38|) we have 

Tf"^^ X^pXre = J c"^^^V^(sAM,y(W)+XA^)X„;3+y" e^"PV^{dx.yA)Xc.f3dy, 

(3.54) 
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where W = w'^^X-re and A = C^^Xro- On the other hand, (W, A) is a solution 
of the foUowing problem 



{-eyiXre^-"" +W,i;) + dy{A,i;))dy = 0, 



(3.55) 



for all (v,V') e Ul,Ay*)m X H^.^iY*) 

Thus (W, A) is a saddle point of the following functional 



defined by 



/(v, i;) = i y c"^^'^\/^(sAM,y (v) + Xx^){so,f},yM + Xap)dy 

+ J" e^°''^^/a{sa0.y{v) + Xa/sjOx^yipdy - ^ J d"^y/ada,y'ipdx.^yi'dy. 

This yields 

/(W,7/i) </(W,A) </(v,A), 

for all (v,V^) e Hi,,(y*)/M x i/i„(r*)/R. 
Consequenthly, for = we get 

/(W, A) > /(W, 0) - ^ y c^/'^^V^ (sa0(v) + X„^)(sA^.,,(v) + > 0. 

Moreover, by taking (v, tp) = (W, A) in (|3.55p we obtain 
c"P^'Xo.f,Xre = 2IiW,A)>0. 



Let us define the function (f> : K"* — > R by 

It is to easy that $ is continuous in endowed with to the norm || ^ || = 
(WW)^- Let 

B e K'^; C symmetric , ||^|| = 1}. 

Since B is compact, $ attains its minimum in B. Then, there exists c > such 
that $ > c in _B, that is 

$ ( ||^°^|| ) > c, for all symmetric tensor ^ = {iaf})- 

This means that c"'^'^^£,af3£,T0 > cS^apirS- The coercivity of (c"^"^^) is now proved, 
(ii) The ellipticity and symmetry of the tensor [d ). 
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First we prove the symmetry. According to (|3.48p we have 

= J [~ey{z'',ya) + dy{y„+'q'^,ya)]dy, (3.56) 

and 

- ey{z^ ,y^)dy = - ey{z^ ,yo. + v'")dy + Cyiz'' ,7j'')dy. (3.57) 
Now, we take (v, V') = (z", ??") in (EHH]). We get 

e^(z^2;„ + 77")dy =/ c,(z^z")dy, 
< -^r (3.58) 

/ e,(z^,7") = / dy{y, + v'',vndy. 

Replacing (|3.57p and p.58p in (|3.56p we obtain 

d"" = j^^ c,(z", Z-) + dy{y^ + ri%y^+ T^^)dy. 

From the above relation, we deduce the symmetry of the tensor (d ). 

Let us prove the coercivity of (d ). Let (X^) be a vector. From (I3.4ip we have 

T'Xa^X^^-j e"^^V^SAM,y(Z)X„dy + y" d'^^V^ {Xx + dxQ)Xo.dy, 

(3.59) 

where Z = tFX^^ and 8 — if X^. It is easy to see that (Z, 8) is the solution of 
the following variational problem 

(c3,(Z,v) +ey(v,X^y<, + 8))dy =0, 
{-ey{7F ,i,)+dy{y„ + jf ,^))dy =0, 

for all (v,'0) e Hi„(y*)/]R x iJi„,(y*)/]R. Moreover, (Z, 8) is a saddle point 
of the following functional 

J : Hi,,(>^*)/M X i/p^,(r*)/R ^ R, 

defined by 



J(v,^) = 




This yields 

J(Z,^) </(Z,8) </(v,8), 
12 




for all (v» G Hi3,(y*)/R x Fi„(F*)/R. 

By taking ^ = in the above inequality, we obtain 

j(z, e) > j(z, 0) = ^ y" (^a + a„^ye)(x, + d„,ye)dy > o. 

On the other hand, by taking (v, tp) = (Z, 6) in (|3.59p we obtain 

T^Xo^X^ = 2/(Z, 9) > 0. 

With the same proof as for the coercivity of {Tf^'^^) we deduce the existence of 
d > such that d XaX^ > dXaXa-- 

The uniqueness of (u, (p) follows by the Lax-Milgram Theorem. 
The proof of Theorem 13. H is now complete. □ 

3.2 Corrector result 

We have the foUwing convergence 

r^(7a/3(u^)) ^7a;3,^(u) + s„/3,y(ui) wcakly in (tj X ^ * ) , (3.60) 

T^iVip^) ^V^ip + Vytp^ weakly in X y*). (3.61) 

The convergence of energies follow easily fro the above relations. Moreover, the 
weak convergences in p.60p and p.6ip are actually strong 

T%jc.0,x{n')) - jc^fiA^) - Sa/3,y(u^) strongly in L\lo x Y*), (3.62) 

T'C^^ip") - V^^ - \7yip^ strongly in L^{uj x Y*). (3.63) 
Now, we can state the following corrector result: 

Theorem 3.2. (correctors). One has the following strong convergences : 

7ai3,x{u^) - lai3,x{u) -W{sai3,y{u^)) strongly in L'^ {uj) , (3.64) 

V ^tp" - V ^tp - W {V yip^) strongly in L^{uj). (3.65) 
Proof. Using convergences p.62p - p.63p and Theorem 12.31 we have 

IcpA^-") - W'(7a/3,4u)) - W'(sa/3,y(ui)) ^ strongly in (3.66) 

V,(p" - Wiy^^) - W{so,pA^^)) ^ strongly in L^{lu). (3.67) 
But 7Q/3,a:(u) G L^(cj) and ^ xV ^ ^^(^) so, by Proposition [3] (i) we get 

W^(7a/3,a;(u)) ^ 7a/3,a;(u) strougly in (tj) , (3.68) 

W[Vxp)^Vx'P strongly in l2(cj). (3.69) 
From (|3.66p - (|3.69p we deduce the convergence p.64p and (I3.65p . □ 
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4 The bending problem 



4.1 The convergence results 

We now consider the variational bending problem 

^1 . \ . ^ f . . ^ (4.70) 

P{xi,X2,z)dz ] v^Va"^ dx + / g'wVa'^dr, 



for all V e W(a;'^), where 



£ e 



W(oj^) = |w G H\lo^)xH'^{uj^)xH'^(uj''); japi^) = and = d^,w^ = on 7° C duj"^ 
We have 

f = (f ) e l2(c.^ X [-1, 1]), q = [q^] e L2(r+ u ri), 



(4.71) 

Due to the ellipticity and symmetry of the bending tensor, by using the Lax- 
Milgram Theorem, we can deduce the existence and uniqueness to solution 
of the variational problem (|4.70p . 
Let us denote ^ 

F\xi,X2) = I f{xi,X2,z)dz. 



J-1 

Thus, the equation (|4.70p takes the form 

2/" c'"/3AM,e^Xc,/3(u")Tp^(v)da; = /" F'v'V^dx+f h'v'V^ dV, 

(4.72) 

for aU V e Wiuj"). 

Theorem 4.1. The sequence (T'^(u^))^ weakly converges to u G W(a;) which 
is the unique solution of the homogenized problem 



2 r ^,a0pa ^^^^^^ ^^^^^^ ^1 1^ r \Y*\a f hV dV , 

•J Juj Jul Jr+ur_ 

(4.73) 

for all V G W(a;), where 



Y* |a= vMy) dy, 
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C""'" = j^^ C^'P-^a {b^^ipe + dl,,y^''^\ dy, (4.74) 
and the local functions w'^^ verifes the local problems 

A {C^PP^^ [S^^S,o + 9^,,,w^«] } ^ ^nu^xY*, ^^^^^^ 
^tB Y* -periodic. 

Remark. We can give another expression for the homogenized bending coeffi- 
cients, that as form : 

C"''"" - j C^'P''^d\y[^'^^ + w"'^] dy, 

where 11"^ = \yay0- 

Proof. To obtain an a priori estimate for u"^, we choose v = u"^ in (|4.72p . It 
follows 

J UJ 

where M > is a positive constant that not depend on e. Using the ellipticity 
of the tensor C^f^P"'" we get 



;/ T^^(u^)dx<A/||u^|U2(, 

Jul" 



Now, the Korn and Poincare's inequalities in perforated domains imply 

I|u1|w(^0 < 

where C does not depend on e. Then, up to a subsequence, (u^) weakly con- 
verges to some u € W((jj). By Theorem 1.2 we deduce that there exists S 
L'^{uj;H^^^(Y*)/R) such that 

T^(u^)^u weakly in L2(tj X r*), (4.76) 
T^Vu') Vu weakly in L^{u; x Y*), (4.77) 

and 

T^V^u") V^u + V^u^ weakly in L^(uj x Y*). (4.78) 
The linearity of implies 

T%T^p{u')) T„;3(u) + dl^^y weakly in L^lo x Y*). (4.79) 

Using now the properties of the unfolding operator , in (|4.72p . we get 

I I T%C"^P^'^V^)T%T^p{u-))T^Tp,i^))dxdy 

= / T^F'v'V^) dxdy + [ T%h'v'V^) dVdy, 
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that is 



= ( T^F'v'v^) dxdy+ [ 



(r5,uri)xy* 



T^h'vWa^) dVdy 



We chose as a test function in (|4.80p 

v^(a;) = vi(x) + £2v2 (x, J) , vi S I?(c^); V2 G 



(4.80) 



Then 



V^v^ {x) = Vvi (x) + e^V^W2 (a^' j) + e^j^va (a;, 



It follows that 



strongly in L^(ijj), 
weakly in L^(ti;), 
V2v^(x) ^ V2vi(x) + V^va (x,?/) weakly in l?luj x F* 



v'^(x) ^ vi(x) 
V,v^(x) ^ Vvi(x) 



And 



r^(Tp,(v^)) - Tp,(vi) + a2^_j^V2 weakly in \.''{uj xY 



(4.82) 



(4.83) 



Using to (I4.82P and ()4.83ll . and passing to the limit in ()4.80|) with e \ 0, we get 



T^p (u) + dlp^yV? Tp„ (vi ) + V2 



dxdy 



= 1 la / F'wi dx+ \Y* \a [ h'v\ dV 



We now let V2 = in (|4.84p . We obtain 



(4.84) 



Ta/3(u) + dlp yU^ Tpa{-vi) dxdy 



Y* U / F'vl dx+ \Y*\a [ h'v\ dV 
Jul Jr+ur_ 



(4.85) 



By density, we deduce 



Y* \a I F'V\ dx+ \Y*\a f h'v\ dV 

Jul Jr+ur_ 



(4.86) 
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for all V G L^(tj; H^^j.{Y*)/'R.). In what follows we chose 

u2 ^ Treiu)^^", w"" e H^^^{Y*)/R). (4.87) 

.2 „,Te 



Then 



Replacing in we obtain 

2 



(4.88) 



= 1 Y* I, / F'v'' dx+ \Y*\a f hV dT, 

(4.89) 



for all V e X2(w; Hl^^{Y*)/R). It follows that 



Tre{u) + T„^(u)92^_^w"'5 Tp,(v) dxdy 



Y* \a I F'V' dx+ \Y*\a hV dr, 

Jr+ur_ 



for all V e L^{lu;H^^^{Y*)/R). This yields 



dy^Tafjin) Tp^(v) dx 



for all V e i2(^.^2^^(-y, 

If we denote 



=1 y* u / F'y' dx+ \Y*\a hV dr, 

Jr+ur_ 

(4.90) 



dy, (4.91) 



by (|4.90p we get the homogenized equation which corresponds to (I4.72p : 

2 f -^Ppa ^^^^^^ Tp„{v)dx -I Y* \a f F'v' dx+ \Y*\a I dT , 

■J Jlj Jul "'r+ur_ 

(4.92) 

for all weL^ito] H^^^{Y*)/R). 

Let us find now the equations verified by the local functions w"^^. 
Letting (j) = in we get 



C<^fip^^ [t„^(u) + dlf,,yV?\ d%^yMx, y) dxdy = 0. 
By density it follows that 

C^fp^iy)^ [t„/3(u) + d^p.yuA dl,,yv{x, y) dxdy = 0, (4.93) 

xY* L J 
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forallveL^(tj;F^„.('^* 
Using now ((4H8)) in we have 



d 

dyp 

and so 

A {C"'3''^^/H [So^r^pe + 5l^,„w-«] } = in X y*. (4.94) 

In order to establish the existence and uniqueness of the solution of (|4.73p . if 
sufRces to prove the coercivity of C in the following sense 

3Ac 7^ Ac(£) > 0, 'i{£,ap)af3 ■ ^af3 = ^fia, C afj-yq^aP^-fri > ■^cS.aP^afS 

Symmetry 

It is easy to check that 



It suffices to prove 



— Q/37JJ _ — /3Q777 _ —afirii 



Starting from the definition (|4.9ip of the C = {C"'^^^), the homogenized bendig 
tensor is evaluated by 

C"'^'"" = J C'^^^V^dlsr.yi^''"' +^^''] dy 

- J^^ C^^'^^ dls^^ylW^ + w^"]a|<.,,,w"'=' dy. (4.95) 
By multiplying (|4.75|) by w"^ and integrating by parts, we prove : 



It follows that 



c 



1^^ C'^'^-'V^ dlsr,y (n^" + ^'''^9)^h,y (n^" + ^''(y)) dy- (4.96) 



From dUHl), we deduce C""''' = C^'"^. 
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EUipticity 

Let {^af3)af3 be a symmetric tensor {^ap — C/3q) ^i^d set 

Using now the coercivity of tensor C'^^"' (x, ?/) and the fact that a ^ 0, we can 
write 

We claim that the second integral in (I4.97P is positive. Assume the contrary. 
Then 

V (5, 77) e {1, 2}2, = ea/j^l,,,^ (n"^ - w"^) - 0. (4.98) 
It follows that 

This implies that 

for some constants and 6, l— 1,2. This yields 

Since ^ ^ 0, we can find an indice (a,/3) such that ^q/j ^ 0. In this case, the 
left-hand side of the above equality is y*-periodic, but the right-hand side is 
not, this is clearly a contradiction. Then the second integral of (|4.97|) is positive 
and so 

C" '^'^ ^ap^'fr, > V i^ap) 7^ Symmetric . 
Let us define the function : R'' R by 

It is easy to see that ^ is continuous in R"* endowed with to the norm 

II T 11= (T„^T„0)i 

Since $ attains its minimum on the unit sphere in R"* and ^' > for all sym- 
metric tensor {^ap) ^ 0, we can conclude that there exists Af > such that 

^'(j^^^) > Af, for all symmetric tensor {iap) ^ 0. 

From the above inequality we deduce 

The uniqueness of the solution of (|4.73p follows now by using the Lax-Milgram 
Theorem. □ 
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4.2 Corrector result 

We have the foUwing convergence : 

T'{Tp,{u')) - Tp,,,(u) + dl^yU2 weakly in h'^iuj x Y*). (4.99) 

The convergence of energies is also proved easily, and implies in particular that 
the weak convergences in (14.99^ is actually strong 

T^(Tp,,,(u^)) ^ Tp,,,(u) + d^p, yU2 strongly in L''{uj x Y*). (4.100) 

Theorem 4.2. (correctors). One has the following strong convergence : 

Tp,,,(u") - Tp„,,{u)-W{d%,yVL2) ^ strongly m l?{uj). 

Proof. We have already seen (see (|4.100|) ) that 

r-'(Tp,,,(u^)) - Tp,,,(u) - d%^yVL2 -> strongly in l.''{uo x F*), 

which, by Theorem l2.3l is equivalent to 

Tp,,,(u^) - W{Tp,^,{n) + d%^yVi2) ^ strongly in l?{uo). 

But Tpo-,a;(u) G L^(w), so from (i) of Proposition [3] on has W{T p^^^i^)) 
TC"p(T,2:(u) strongly in L^(tj), whence the desired result. □ 

5 Conclusion 

In this paper we have rigously established the limiting equations modelling the 
behavior of a thin piezoelectric perforated shells, i.e., we have explicity described 
the homogenized coefficients of the elastic, dielectric and coupling tensors (for 
details, see [9]). 
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